THE  5th  ASIAN  COMPUTAITIONAL  FLUID  DYNAMICS 
BUSAN,  KOREA,  OCTOBER  27-30, 2003 


DNS  of  turbulent  boundary  layer  with  time-periodic  blowing 

through  a  spanwise  slot 

K.  Kim  1  and  H.J.  Sung  2 

1.  Department  of  Mechanical  Engineering ,  Korea  Advanced  Institute  of  Science  and  Technology , 

373-1,  Guseong-dong,  Yusong-gu,  Daejeon,  305-701,  Korea, 
e-mail:  kyoungyoun@webmail.kaist.ac.kr 

2.  Department  of  Mechanical  Engineering,  Korea  Advanced  Institute  of  Science  and  Technology, 

373-1,  Guseong-dong,  Yusong-gu,  Daejeon,  305-701,  Korea, 
e-mail:  hjsung@kaist.ac.kr. 

Corresponding  author  H.J.  Sung 


Abstract 

Direct  numerical  simulations  are  extended  to  see  the  effects  of  time-periodic  blowing  frequency 
from  a  spanwise  slot  on  a  turbulent  boundary  layer.  The  time -periodic  blowing  frequency  is  given  in  a 
range  0<f  <0.08  at  a  fixed  blowing  amplitude  of  A+= 0.5.  The  effects  of  the  blowing  frequency  are 
scrutinized  by  examining  the  phase-  or  time-averaged  turbulent  statistics.  A  most  effective  blowing 
frequency  is  obtained  at  f= 0.03,  where  the  maximum  increases  of  Reynolds  shear  stress,  streamwise 
vorticity  fluctuations  and  energy  redistribution  are  made.  The  phase-averaged  stretching  and  tilting 
terms  are  studied  for  analyzing  the  increase  of  streamwise  vorticity  fluctuations  which  is  closely 
related  to  turbulent  coherent  structures. 

Keyword:  periodic  blowing,  spanwise  slot,  turbulent  boundary  layer,  DNS,  energy  redistribution 

1.  Introduction 

Advances  in  the  understanding  the  coherent  structure  of  wall-bounded  turbulent  flow  have 
intensified  interest  in  controlling  the  near-wall  turbulence.  Recent  direct  numerical  simulations  have 
demonstrated  a  successful  active  control  over  the  entire  wall  [1,2].  From  a  practical  point  of  view, 
however,  this  active  control  over  the  entire  wall  is  difficult  to  implement  because  it  requires  a  dense 
population  of  sensors  and  actuators  on  the  wall.  Many  attempts  have  been  made  to  devise  a  practical 
method  for  controlling  wall  bounded  flows.  These  include  the  modification  of  the  wall  surface  by 
installing  riblets  [3],  as  well  as  the  use  of  a  compliant  wall  [4],  a  wall  deformation  [5]  or  a  spanwise 
oscillating  wall  [6].  Among  the  approaches  considered  to  date,  the  use  of  local  suction/blowing  [7,8] 
deserves  more  detailed  study  because  it  provides  an  efficient  and  simple  means  for  locally  actuating 
the  wall-bounded  flow.  Moreover,  the  strength  of  the  actuation  can  be  controlled  with  relative  ease  by 
local  suction/blowing. 

Most  previous  experimental  and  numerical  studies  of  local  suction/blowing  have  focused  on 
steady  actuation  [9,10,11,12].  It  is  reported  that  the  local  steady  blowing  lifts  up  near-wall  streamwise 
vortices,  thereby  reducing  the  interaction  of  the  vortices  with  the  wall.  The  steady  blowing  leads  to  a 
reduction  in  the  skin  friction  near  the  wall,  combined  with  an  increase  in  the  turbulent  intensity  and 
skin  friction  far  downstream  from  the  slot.  In  contrast  to  the  previous  studies  that  considered  only 
steady  blowing,  a  relatively  few  studies  of  unsteady  suction/blowing  were  made  experimentally  and 
numerically  [13,14,15,16,17].  Park  et  al.  [13]  performed  experiments  to  probe  the  effects  of  periodic 
blowing  and  suction  through  a  spanwise  slot  on  a  turbulent  boundary  layer.  The  higher  forcing 
frequency  induces  greater  changes  in  the  turbulent  structures  of  boundary  layer.  Rhee  and  Sung  [14] 
performed  unsteady  Reynolds-Averaged-Navier-Stokes  simulations  and  compared  the  simulation 
results  with  those  of  the  experiments  of  Park  et  al  [13].  It  is  known  that  the  near- wall  streamwise 
vortices  play  a  dominant  role  on  the  wall-bounded  flows  [18],  however,  the  responses  of  near- wall 
vortices  to  the  unsteady  periodic  blowing  were  not  studied  in  detail. 
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Fig.  1.  Schematic  diagram  of  the  computational  domain 

Recently,  Kim  and  Sung  [19]  performed  direct  numerical  simulations  to  analyze  the  effects  of 
time -periodical  blowing  through  a  spanwise  slot  on  a  turbulent  boundary  layer.  The  blowing  velocity 
was  varied  in  a  cyclic  manner  from  0  to  2 A+  (A+= 0.25,  0.50  and  1.00)  at  a  fixed  blowing  frequency  of 
/=0.017.  The  effect  of  steady  blowing  (^4+=0)  was  also  examined,  and  the  steady  blowing  results  were 
compared  with  those  of  unsteady  periodic  blowing.  The  periodic  blowing  generates  a  spanwise 
vortical  structure  in  the  downstream  of  the  slot.  The  energy  redistribution  is  enhanced  by  the  periodic 
blowing.  However,  it  is  expected  that  the  turbulence  structure  is  more  sensitive  to  the  blowing 
frequency  than  the  blowing  strength. 

In  the  present  study,  the  effect  of  the  blowing  frequency  on  turbulent  boundary  layer  is  studied  as 
a  sequel  to  the  previous  study  of  Kim  and  Sung  [19].  Main  emphasis  of  this  study  is  placed  on  the 
blowing  frequency  effect  on  near-wall  turbulent  flow  structures  downstream  of  the  spanwise  slot.  The 
Reynolds  number  based  on  the  momentum  thickness  at  inlet  is  Reg= 300,  and  the  slot  width  is 
approximately  100  wall  units.  The  localized  time -periodic  blowing  is  given  by  changing  the  vertical 
velocity  on  the  spanwsie  slot.  The  blowing  frequency  is  in  a  range  of  0</^ <0.08  at  a  fixed  blowing 
amplitude  (^4+=0.5).  The  frequency  responses  are  scrutinized  by  examining  the  phase-  or  time- 
averaged  turbulent  statistics.  A  most  effective  blowing  frequency  is  observed  at  7^=0. 03,  where  the 
maximum  increases  of  Reynolds  shear  stress,  streamwise  vorticity  fluctuations  and  energy 
redistribution  are  made. 

2.  Computational  details 

In  the  present  study,  the  computational  details  are  the  same  as  those  of  the  earlier  work  of  Kim 
and  Sung  [19].  The  flow  configuration,  boundary  conditions  and  other  numerical  procedures  are 
summarized  in  the  following.  As  shown  in  Fig.  1,  the  domain  size  is  2006^x30 6^x40 6in  in  the 
streamwise,  wall-normal  and  spanwise  directions,  where  the  corresponding  mesh  size  is  257x65x129. 
The  mesh  is  uniform  in  the  streamwise  and  spanwise  directions,  but  a  hyperbolic  tangent  stretching  is 
used  in  the  normal  direction  to  cluster  points  near  the  wall.  The  mesh  resolutions  are  Zlx+«12.40, 
Ay+min~<).\l ,  Ay+max& 23.86,  and  Zlz+»4.96,  based  on  the  friction  velocity  at  the  inlet.  Realistic  velocity 
fluctuations  at  the  inlet  are  obtained  using  the  method  of  Lund  et  al.  [20].  The  convective  outflow 
condition  (duildt)+c(duildx)  =  0  is  used  at  the  exit,  where  c  is  taken  to  be  the  mean  exit  velocity.  A  no¬ 
slip  boundary  condition  is  imposed  at  the  solid  wall.  At  the  free-stream,  the  conditions  u=Uoo  and 
dv/dy=dw/dy= 0  are  imposed.  Periodic  boundary  conditions  are  used  in  the  spanwise  direction.  The 
spanwise  slot  for  periodic  blowing  extents  from  x=75.80in  to  x=82.0  6in,  where  the  location  of  the  inlet 
is  defined  as  x=0.  The  slot  width  is  6+«100  in  wall  units.  The  periodic  blowing  at  the  slot  is  generated 
by  varying  the  wall-normal  velocity  according  to  the  equation: 

^  slot  !U«,=A{\  +  COS2  7ft)  (1) 
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The  maximum  blowing  velocity  (vsiot=2 A)  is  imparted  at  t=0/4T  and  the  minimum  (v^=0)  at  t=2/4T ' 
where  T  is  the  blowing  period.  At  t=M4T  and  3/4 T,  the  blowing  velocities  are  the  same  as  that  of 
steady  blowing  with  decelerating  and  accelerating  phase,  respectively.  The  amplitude  of  periodic 
blowing  is  A+= 0.5  in  wall  unit,  which  corresponds  to  the  value  of  vrms  aty+=15  without  blowing.  The 
blowing  frequency  (f=fv/uT>in2)  varies  in  a  range  0 <f  <0.08,  where  uTjn  is  the  friction  velocity  at  the 
inlet.  The  governing  Navier-Stokes  and  continuity  equations  are  integrated  in  time  by  using  a 
fractional  step  method  with  an  implicit  velocity  decoupling  procedure  [21].  A  second-order  central 
difference  scheme  is  used  in  space  with  a  staggered  mesh.  The  Reynolds  number  based  on  the 
momentum  thickness  at  the  inlet  is  Reg= 300.  The  computation  time  step  is  AtUJ 6{„=0.3,  which 
corresponds  to  At+& 0.25  in  wall  units.  The  total  time  over  which  statistical  averages  are  calculated  is 
Tavg=\%0006in/Uoo,  which  corresponds  to  150,  460  and  1250  periods  for  f= 0.01,  0.03  and  0.08, 
respectively.  The  imposition  of  periodic  blowing  may  lead  to  periodic  variations  in  the  global 
physical  quantities  of  the  flow.  Hence,  it  is  necessary  to  represent  each  flow  quantity  as  a 
superposition  of  three  components 

q(x,  y,z,t)  =q(x,y)  +  q  (x,  y,t)  +  q"(x,y,z,t),  (2) 

where  the  instantaneous  quantity  q  is  decomposed  into  a  time-mean  component  q  ,  an  oscillating 
component  q  and  a  random  fluctuating  component  q".  The  time-average  is 


q(x,y)  = 


- —  w  q(x,  y,  z,  t)dzdt , 


TtotLz 


(3) 


where  Ttot=NT  is  the  time  over  which  the  quantity  is  averaged  and  N  is  the  total  number  of  periods. 
The  oscillating  component  q  is  obtained  from  the  relation 

q(x,y,t)=<q>(x,y,t)-q(x,y )  (4) 

where  <q>(x,y,t )  is  the  phase-average,  which  is  defined  as 


Accordingly,  the  random  fluctuation  component  q"  is  expressed  as 

q"(x,y,z,t)  =  q(x,y,z,t)-<q>  (x,  y,  t )  (6) 


<  q  >  (x,y)  = 


NL 


i  fl, 

2J0  q(x,y,z,t)dz. 

n= 1 


(5) 


3.  Results  and  Discussion 

It  is  important  to  find  the  energy  redistribution  by  the  present  unsteady  blowing.  The  relative 
increase  of  turbulent  intensity  is  defined  as  Aq(%)  =  {qmax-q max,o)  /qmax,ox  100,  where  the  subscript  "o" 
denotes  no  blowing.  As  displayed  in  Fig.  2,  the  increases  of  cross-stream  components  Av"nns  and 
Aw"rms  are  larger  than  that  of  streamwise  component  Au"rms.  This  means  that  the  intercomponent 
energy  transfer  is  enhanced  by  unsteady  blowing.  The  behaviours  of  Av"nns  and  Aw"rms  indicate  that 
they  have  a  broad  peak  at  around /*=0.03.  Next,  to  see  the  relative  contribution  to  the  turbulent  kinetic 
energy  of  streamwise  turbulent  intensity  and  intensities  normal  to  the  mean  flow,  the  energy  partition 
parameter  K*  is  calculated,  K*=2u"2 l(y"2+w"2)  [22].  The  parameter  K*  is  measured  at  which  the 
increase  of  turbulent  kinetic  energy  is  maximum.  It  is  obvious  that  K*  has  a  local  minimum  at  f= 0.03. 
The  smaller  value  of  K  means  that  the  energy  redistribution  is  more  active.  However,  as  the  blowing 
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Fig.  2.  Relative  increases  of  turbulent  intensities  and  the  energy  partition  parameter  JC 


Fig.  3.  Time-averaged  values  of  (a)  skin  friction  coefficient,  (b)  wall  pressure, 

(c)  wall  pressure  fluctuations. 

frequency  increases  further  (/^=0.08),  K*  converges  to  that  of  '  f= 0.  This  suggests  that  the  flow  is  not 
sensitive  to  the  higher  frequency  blowing. 

Figure  3  shows  the  streamwise  distributions  of  time-averaged  skin  friction  Cj,  wall  pressure  pw 

and  rms  of  wall  pressure  fluctuations  p"w,rms •  Three  cases  of  the  blowing  frequency  (/*=0.01,  0.03  and 
0.08)  are  chosen  to  distinguish  the  flow  structures  by  unsteady  blowing.  The  effective  blowing 
frequency  is  chosen  at /*=0.03,  which  gives  the  minimum  value  of  K*.  The  lower  and  higher  blowing 
frequencies  are  chosen  at /*=0.01  and /*=0.08,  respectively.  /*=0  corresponds  to  'steady  blowing'  and 
A+= 0  is  'no  blowing'.  For  all  blowing  cases  in  Fig.  3(a),  rydecreases  rapidly  near  the  slot  and  increases 
in  the  downstream.  It  is  seen  that  the  recovery  of  Cf  at /^=0.03  is  faster  than  other  two  blowing  cases. 


838 


DNS  of  turbulent  boundary  layer  with  time-periodic  blowing  through  a  spanwise  slot 


Fig.  4.  Time-averaged  streamwise  velocity  profiles 


f*=Q.Ql  f*=0M3  f+=0.08 
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x/0m  x/8j„  l/ftn. 

Fig.  5.  Contours  of  <cq7>-cdzo  .  The  contour  levels  are  from  -0.25  to  0.25  by  increments 

However,  the  variations  of  wall  pressure  pw  in  Fig.  3(b)  are  almost  the  same  regardless  of  the 

blowing  frequency.  The  adverse  pressure  gradient  appears  ahead  and  behind  the  slot,  whereas  the 
favourable  pressure  gradient  occurs  above  the  slot.  As  shown  in  Fig.  3(c),  the  wall  pressure 
fluctuations  p"w>rms  are  very  sensitive  to  the  effective  blowing  frequency  (/*=().  03). 

Distributions  of  the  time-averaged  streamwise  velocity  at  four  locations  (x/Oin= 80,  90,  100  and 
110)  are  displayed  in  Fig.  4.  As  compared  to  the  case  of  'no  blowing',  a  region  of  retarded  flow  is 
observed  near  the  wall.  As  the  flow  moves  downstream,  the  region  of  retarded  flow  gradually  shifts 
away  from  the  wall  and  finally  decays.  However,  the  time-averaged  streamwise  velocity  is  invariant 
with  the  blowing  frequency.  This  is  consistent  with  the  previous  results  [16,17]. 

Contours  of  the  difference  between  the  phase-averaged  spanwise  vorticity  for  periodic  blowing 
<coz>  and  the  time-averaged  spanwise  vorticity  for  'no  blowing'  ~coz  0  are  shown  in  Fig.  5  during  one 

period  (17).  A  negative  region  of  <coz>-(Dz  0 ,  denoted  by  dashed  line,  appears  above  the  slot  for  three 
cases.  This  is  because  a  negative  spanwise  vorticity  layer  in  the  vicinity  of  the  wall  is  shifted  upward 
by  blowing.  Note  that  the  negative  (positive)  value  of  <coz>-coz  0  represents  the  increase  (decrease)  of 

the  magnitude  of  spanwise  vorticity  \coz\.  This  is  caused  by  the  fact  that  cozo  is  negative  inside  the 

boundary  layer.  For  /^O.Ol,  it  is  clearly  seen  that  a  region  of  strong  negative  spanwise  vorticity  is 
formed  above  the  slot  and  convects  downstream  as  time  goes  by.  During  the  accelerating  phase 
(*=3/47M)/47),  the  formation  of  the  strong  negative  spanwise  vorticity  occurs  above  the  slot,  which 
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f+=0.01 


f+=om 


f+=om 


Fig.  6.  Contours  of  the  difference  between  the  time-averaged  turbulent  intensities  and 
Reynolds  shear  stress  for  periodic  blowing  and  those  for  no  blowing. 


results  from  the  lifted  wall  vorticity  layer  by  blowing.  The  region  of  strong  negative  spanwise  vorticity 
convects  to  downstream  during  the  decelerating  phase  (t=\/4T~2/4T),  when  the  adverse  pressure 
gradient  decreases  above  the  slot.  For  f= 0.03,  the  lifted  vorticity  layer  convects  to  lesser  extent  as 
compared  with  /*=0.01.  This  may  be  attributed  to  the  fact  that  the  time  difference  between  the 
accelerating  and  decelerating  phases  becomes  smaller  with  increasing  the  blowing  frequency.  Thus,  a 
newly  generated  strong  spanwise  vorticity  coexists  with  the  weaker  prior  one  which  convected  in  the 
decelerating  phase  of  the  previous  period.  For  /*=0.08,  however,  the  time  difference  is  so  small  that 
the  afore-stated  unsteady  responses,  such  as  the  convection  of  the  strong  negative  spanwise  vorticity, 
are  not  found. 

Figure  6  shows  contours  of  the  difference  between  the  time-averaged  rms  velocity  fluctuations 
(u" rms ,  V'rms,  w" rms  and  -w"v")  for  blowing  and  the  values  of 'no  blowing'.  For  three  blowing  cases,  the 
turbulent  intensities  and  Reynolds  shear  stress  are  significantly  enhanced  downstream  of  the  slot.  The 
maximum  increase  of  u”rms  is  located  closer  to  the  slot  than  those  of  v"rms,  w"rms  and  -w"v".  The  general 
contour  patterns  for  three  blowing  cases  look  similar.  However,  a  closer  inspection  of  the  patterns 
reveals  that  the  maximum  increases  of  v"rms,  w"rms  and  -u”v"  for  f= 0.03  are  located  closer  to  the  slot 
as  compared  with  other  blowing  frequencies.  This  may  be  related  to  the  rapid  increase  of  c/behind  the 
slot  for /=0.03,  as  mentioned  earlier  in  Fig.3(a). 

Profiles  of  the  time-averaged  rms  velocity  fluctuations  and  Reynolds  shear  stress  are  shown  in 
Fig.  7.  The  velocity  fluctuations  and  Reynolds  shear  stress  for  blowing  are  increased  downstream  of 
the  slot  compared  with  the  case  of  'no  blowing'.  The  maximum  increase  of  u  "rms  is  located  at  around 
x/0in=  100.  However,  the  maximum  increases  of  v”rms,  w"rms  and  -u”v"  are  located  in  the  farther 
downstream  (x/0in=  110).  It  is  interesting  to  find  that  the  increase  of  u"rms  for /*=(). 03  is  smaller  that 
those  for  other  frequencies,  whereas  the  increases  of  v"rms,  w”rms  and  -u”v"  for  0.03  are  larger  than 
those  for  other  frequencies.  This  suggests  that  /*=0.03  is  the  most  effective  blowing  frequency  in 
promoting  the  energy  redistribution.  Park  et  al.  [13]  reported  that  the  transverse  fluctuating 
components  are  more  susceptible  to  disturbances  than  the  streamwise  component.  On  the  other  hand, 
the  time-averaged  turbulent  intensities  appear  to  be  insensitive  to  the  higher  blowing  frequency 
(7^=0. 08). 

Figure  8  shows  the  time-averaged  rms  vorticity  fluctuations  (a>i"rms).  All  components  of  the 
vorticity  fluctuations  (cox"rms,  coy"rms  and  coz"rms)  are  enhanced  downstream  of  the  slot  by  periodic 
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Fig.  7.  Time-averaged  turbulent  intensities  Fig.  8.  Root  time-mean  square  of  vorticity 
and  Reynolds  shear  stress.  fluctuations. 

blowing.  A  closer  examination  of  Fig.  8  indicates  that  coy"rms  and  coz"rms  are  enhanced  by  blowing 
regardless  of  the  blowing  frequency.  However,  coxnnns  is  most  enhanced  at  a  particular  blowing 
frequency  (Z^O.03).  Accordingly,  the  most  effective  blowing  frequency  is  /*=0.03,  at  which  the  near¬ 
wall  vortical  structure  is  most  activated.  It  is  known  that  the  location  of  local  maximum  of  cox"rms 
corresponds  to  that  of  the  center  of  the  streamwise  vortex  in  the  wall  region  [23].  Recall  that  the 
maximum  increases  of  wall  pressure  fluctuations  and  Reynolds  shear  stress  are  observed  at  /^=0.03 
(Figs.  3  and  7). 

Conclusion 

Detailed  numerical  analysis  has  been  performed  to  see  the  effect  of  blowing  frequency  on  a 
turbulent  boundary  layer.  The  slot  width  is  Z?+»100  in  wall  units  and  the  blowing  frequency  varies  in  a 
range  0<f  <0.08  with  a  fixed  blowing  amplitude  (^4+=0.5).  An  effective  blowing  frequency  is  observed 
at /=0.03,  which  gives  the  minimum  value  of  K*.  The  time-averaged  streamwise  velocity  and  wall 
pressure  are  invariant  with  the  blowing  frequency.  However,  the  time-averaged  skin  friction  and  rms 
of  wall  pressure  fluctuations  are  sensitive  to  the  blowing  frequency.  Furthermore,  the  recovery  of  skin 
friction  is  fast  at  f= 0.03.  The  maximum  increase  of  u "rms  is  located  closer  to  the  slot  than  those  of 
v'U,  W'rms  and  -w"v".  The  increase  of  u ”rms  for  /^=0.03  is  smaller  than  those  for  other  frequencies, 
whereas  the  increases  of  v"rms,  w"rms  and  -u"v"  for  f= 0.03  are  larger  that  those  for  other  frequencies. 
The  time-averaged  streamwise  vorticity  fluctuations  cox"rms  is  most  enhanced  at  ^=0.03,  which 
activates  the  near-wall  vortical  structures. 
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